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Abstract. We present the results of combined optical and near-infrared photometry for the globular cluster systems 
l/~) ' of the giant ellipticals IC 4051 and NGC 3311. We use the reduced age-metallicity degeneracy in (V — I) vs. (V — H) 

color-color diagrams to derive the cumulative age distribution within the red sub-population of globular clusters 
and to search for age sub-populations. The age distribution is then compared to the one determined for simulated 
globular cluster systems in order to set constraints on the relative age and size of these globular cluster sub- 
populations. In both galaxies we find a significant fraction of globular clusters with ages between 2- 5 Gyr. We 
O !' also investigate the metallicity distribution in both systems. Small number statistics prevent us from making any 

definite statements concerning NGC 3311, but we find that the derived metallicity distribution of the IC 4051 
clusters strongly depends on the assumed age distribution. Based on our most likely result that finds a large number 
of young/intermediate age clusters (^2 Gyr) within the selected globular cluster sample, we find metallicity peaks 
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at ~ -0.2 for the old clusters and +0.8 for the young clusters. Only few very metal poor clusters are found. However, 
the metallicity distribution within the young/intermediate globular cluster population is significantly affected by 
our choice of the applied Single Stellar Population model. The mean metallicity of the second generation of 
globular clusters changes from the above mentioned and extremely high +0.8 dex to +0.2 dex. Note that the 
■ model dependency becomes less severe with an increasing age of the cluster population. 
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Key words, globular clusters: general — globular clusters: individual (NGC 3311, IC 4051), sub-populations, age 
distribution, Monte-Carlo simulations 



1. Introduction 

Globular cluster systems (GCSs) are widely used 
to probe the formation and evolution of their host 
galaxies (e.g. ISearle fc Zinn 19781 JZepf fe Ashman 1993| 
|Ashman fc Zepf 1998| |Kissler-Patig et al. 2002| ) although 
they represent only a small fraction of the galaxy's to- 
tal luminosity. Hosting an almost perfect single stel- 
lar population (SSP), globular clusters allow the in- 
vestigation of galaxies which, due to their large dis- 
tances, cannot be resolved into single stars. The inves- 
tigation of massive early-type galaxies, not found in the 
neighborhood of the Milky Way, in particular suffers 
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from the large distance drawback. The diffuse light of 
a galaxy is a composite of all underlying stellar pop- 
ulations which can not be disentangled, and popula- 
tions of several ages and/or metallicities might indeed be 
hidden in the integrated light (e.g. ILarsen et al. 20 03). 
Detailed information on the age structure of the galaxy 
is therefore very difficult to obtain from the diffuse in- 
tegrated light. Globular clusters, on the other hand, are 
luminous enough to be observed even at distances of 
~100 Mpc. Their stars share the same age and metallic- 
ity, and especially early-type galaxies can host quite nu- 
merous GCSs ( |Ashman fc Zepf 1 998|). The availability of 
SSP models (e ' g. IBruzual 20001 IBruzual fc Chariot 20031 
IVazdekis 19991 IVazdekis et al. 20031 IMaraston 2005|) for 
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a wide range of ages and metallicities, applied to deep 
photometric or spectroscopic data of GCSs, allows us 
to set constraints on their age and chemical com- 
position (e.g. [Uoudfrooij et al. 2001a| IPuzia et al. 2 002 
IPuzia et al. 2004bl |Hempel et al. 2003| > . 

This basic tool allows extragalactic GC research 
to probe a variety of important questions. Since the 
discovery of the bimodal color distributions in globu- 
lar cluster systems, their origin has been a source of 
ongoing debate, although there is general agreement 
that distinct globular cluster populations of different 
age and/or metallicity are the origin of this bimodality. 
When and how these populations were formed and what 
can be inferred for the evolution of the host galaxy 
are some of the most outstanding questions addressed 
today, with crucial implications for galaxy formation 
theories. Consequently, much effort has been put into 
the observation of GCSs in early- type galaxies, which 
were classically assumed to have formed early on and 
to have evolved passively since then. From previous 
surveys, it is known that early-type galaxies are, on the 
one hand, well represented by the fundamental plane, 
indicating basic similarities, whereas their GCSs can 
show significant differences; for example with regards 
to metalli city (e.g. [G eisler, Lee fe Ki m 1996| ), color dis- 
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Har ris 2001(1 . Therefore, much caution has to be paid 
before claiming a universal formation scenario for early- 
type galaxies. In fact, it seems more and more likely 
that we have to abandon the idea of a single formation 
theory and find external parameters which might be 
decisive for the particular formation and evolution of a 
given galaxy. Factors involving the galaxy environment, 
such as its position within the host galaxy cluster/group 
or its isolation from other galaxies, is likely to be one 
of the most important external parameters. Galaxy 
interactions (e.g. mergers and accretion events) seem to 
be an important, if not the dominant, process, governing 
galaxy formation and evolution (see contributions in 
IWielen 19 90). Consequently, the GCSs hosted by galaxies 
which were likely involved in such interactions need to be 
investigated. First results on the investigation of early- 
type galaxies in low-density environments were published 
by Kuntschner et al. (2002), who used integrated- light 
spectroscopy to derive ages and metallicities of 40 E and 
SO galaxies and conclude that early-type galaxies in a 
low density environment show a broad distribution of 
luminosity-weighted ages and are on average younger 
than cluster ellipticals. Unfortunately, no early-type 
galaxies in a high galaxy density environment exist in 
the very local Universe. The determination of globular 
cluster ages by integrated-light spectroscopy is therefore 
quite challenging due to the faintness of the targets. 
The capability of high signal-to-noise spectroscopy with 
respect to age and metallicity determination has been 
demonstrated by various studies (e.g. IPuzia et al. 20 04a 
IPuzia et al. 2004bj> and spectroscopy will doubtless be 



the method of choice when precise age estimates of a 
subsample of bright clusters are required. However, as 
long as the objective is to detect different broad age 
subpopulations in GCSs and to set coarse age limits 
on the GCS as a whole, then photometric studies (e.g. 
IKundu fc Whitmore 2001al IRuiidu fc Whitmore 2001 bl 
IGebhardt fe Kissler-PatTg~1 999) are the most efficient 
alternative. 

Until recently, the major drawback for the use 
of photometric studies to investigate the age distri- 
bution in GCSs was the age-metallicity degeneracy 
QWorthey~ 994). With the introduction of combined opti- 
cal and near-infrared photometry fe.g. lMinniti et al~ 1996 
IGoudfrooij et al. 2001b| |Kissler-Patig et al. 20021 it is 
now possible to at least partially overcome this obsta- 
cle and derive the relative ages of globular cluster sub- 
populations, as well as their metallicities. The extended 
wavelength range allows one to recognize intermediate age 
globular clusters (e.g. those formed in a recent merger) 
from the bulk of first generation objects, formed shortly 
after the Big Bang. In addition, colors like V-H provide 
a sensitive metallicity index. The integrated IR colors for 
Galactic globular clusters of Cohen & Matthews (1994) 
show that the V — H index has three times the metallicity 
sensitivity of V — I, making it a good index for investigat- 
ing the reality of any multiple metallicity populations. 

In this study we investigate the GCSs of two important 
early-type galaxies, viz. NGC 3311, representing a giant 
cD galaxy at the center of a high density environment (the 
Hydra I galaxy cluster), and IC 4051, a giant E2 in the 
Coma galaxy cluster with no luminous neighbor, located 
several cluster core radii from the center. This will pro- 
vide a range of environmental densities in order to begin 
to probe their effect on the GCS. The GCS of NGC 3311 
has been investigated in the optical by many different 
groups (e.g. lGrillmair et al. 199 4 McLaughl in et al. 1 995 
|Brodie, Larsen, fc Kissler-Patig 2000| ISecker et al. 199511 . 
They found a populous GCS which appears to show 
the general color bimodality of such systems. Due 
to its central position in a high density environ- 
ment, several wide field surveys have been car- 
ried out (e.g. |Harris, Smith, fe Myra 1983| IHTlker 2003al 
Hilker 2003b) to investigate spatial features in the 
GCS. NGC 3311 is of special interest in globular 
cluster studies, since it is not only the third clos- 
est cD galaxy, following in distance M 87 in Virgo 
and NGC 1399 in Fornax, but also hosts one of the 
richest GCSs known (e.g. |Harris, Smith, fc Myra 1 983 
|McLaughlin et al. 1995| ). The GCS of IC 4051 has only 
been attainable until now from HST imaging and was the 
subject of papers by Baum et al. (1997) and Woodworth 
and Harris (2000). The latter work used V, I photometry 
and found tentative evidence of bimodality, with metal- 
intermediate ([Fe/H]~ —0.7) and metal-rich (~ solar) 
peaks and a near-complete lack of metal-poor clusters. 
The lack of metal-poor clusters and the very high specific 
frequency (Sn = 11±2) that they found makes the GCS of 
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Property 


NGC 3311 


IC 4051 


Reference 


"r> a / TonrirA 
KA(J2000) 


lUh 3om 43s 


Ion UUm 55s 


deVaucouleur et al. 19911 


DEC(J2000) 


-27° 31' 42" 


+28° 00' 27" 


deVaucouleur et al. 199ll 


Bt,o 


12.65 


14.17 


de Vaucouleur et al. 19911 


Eb-v 


0.079 


0.011 


Schlcgcl, Finkbeiner & Davis 1998 




A v 


0.263 


0.035 


Schlegel, Finkbeiner & Davis 1998 




A { 


0.154 


0.02 


Schlegel, Finkbeiner & Davis 1998 




A H 


0.046 


0.006 


Schlegel, Finkbeiner & Davis 1998 




(B - V)e f/ 


1.04 


1.01 


deVaucouleur et al. 19911 


(m - M) v 
M v 


33.70 ± 0.08 
-22.07 ± 0.20 


34.82±0.13 
-21.62 ± 0.35 


Jensen, Tonry & Thomson (2001)| 
deVaucouleur et al. 199ll 



IC 4051 very unique and worthy of further study. General 
information concerning these 2 galaxies can be found in 
Table 1. 

In this paper, we combine optical and near-infrared 
photometry to derive the cumulative age distribution (e.g. 
|Hempel k Kissler-Patig 2004| | which we will compare to 
simulated GCSs of known age structure (size and age 
of sub-populations). In this particular case the model 
isochrones by Bruzual & Chariot (2003) have been applied 
and the results are discussed in great detail. Nevertheless, 
an additional section will be spent to discuss the depen- 
dency of the results on the choice of SSP models. To do so 
we will repeat the analysis and apply the latest version of 
the SSP models by Maraston (2005). These results enable 
us to set constraints on the age structure in both systems 
and to infer clues to the evolutionary history of their host 
galaxies ( |Hempel fc Kissler-Pati g 2004 ) . The paper is or- 
ganized as follows. In Section |21 we describe the observa- 
tions as well as the data reduction. The derivation of the 
cumulative age distribution within the observed samples is 
described briefly in Sect. 13.21 For a detailed description we 
refer to Hempel et al. (2003) and Hempel & Kissler-Patig 
(2004). In Section 0] we derive and discuss the metallic- 
ity distribution. We summarize our results in Section |3] 
We will dedicate Appendix ^] to the discussion of how 
the results are influenced by our choice of a Single Stellar 
Population model. 

2. Observations and Data Reduction 

The optical and the near-infrared data for both galaxies 
were reduced by different groups, using different software. 
Therefore we will give a brief description for each data set. 

2.1. Optical data: HST/WFPC2 

Both targets were observed with HST + WFPC2 in two 
optical pass bands, NGC 3311 in F814W and F555W, and 
IC 4051 in F814W and F606W, respectively. NGC 3311 
observations are part of a large survey on extragalac- 
tic GCSs, carried out under program number GO. 6554. 
The observations of IC 4051 were carried out under pro- 



Table 2. Exposure times for NGC 3311 and IC 4051 



Filter (Instrument) 


NGC 3311 


IC 4051 


F814W (WFPC2) 


3800 sec 


5200 sec 


F555W (WFPC2) 


3700 sec 




F606W (WFPC2) 




20500 sec 


F160W (NICMOS2) 


2560 sec 


2560 sec 



gram GO. 6283, dedicated to stellar populations in ellipti- 
cal galaxies. The total exposure times in the various pass 
bands are given in Table El for both targets. 
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Fig. 1. Comparison between (V — I) colors based on aper- 
ture and psf photometry. The solid line represents the 
least-squares fit to (V — I) based on aperture photom- 
etry. The objects in this plot contain the (V — I) colors 
for the PCI chip as well as for all three wide field chips. 
As mentioned in Sect. El the latter could not be included 
in our age analysis, due to the smaller spatial coverage of 
the near- infrared observations. 

For IC 4051 we took the photometry directly from 
Woodworth & Harris (2000). The reduction of the optical 
data for NGC 3311 is described in the following. 
The photometry on NGC 3311 clusters was done us- 
ing the DAOPHOT package for aperture- and alterna- 
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tively for psf photometry as well as using SExtractor 
(Bcr tin fe Arnout 1996(1 aperture photometry. In Figure ^ 
we show the comparison between (V — I), calculated from 
instrumental V and I band magnitudes, using DAOPHOT 
p,s/-photometry and SExtractor aperture-photometry for 
all objects detected on the PCI chip, for which later on 
also H-band data were obtained. The extremely red ob- 
jects ((V — I)> 1.6) have been rejected for being objects 
either in the vicinity of the central dust shell found in the 
galaxy core or close to the edge of the PC chip. Including 
an error selection, i.e. S(V — I)<0.15 mag, the results of 
a least square fit between the (V — /) results are: 

(V - I)aperture = 0.983(F - I) ps f + 0.015 (1) 

There is no difference between the results, the color 
offset being smaller than the photometric error, except 
for very red objects ((V - /)> 2.0 mag). The NGC 3311 
results presented in this Paper are based on SExtractor 
photometry. 

The instrumental magnitudes were transformed to 
the Johnson V and Cousin I magnitudes according to the 
procedure given by Holtzman et al. (1995) and reddening 
corrected using the values given in Tabled The aperture 
corrections for both filter bands were determined as the 
mean difference (21 objects) between a 6 pixel instrumen- 
tal magnitude and the 075 aperture used by Holtzman 
et al. (1995). The calculated corrections, 0.229 mag and 
0.265 mag for the V and I band, respectively, were than 
applied to all of the V- and I-band detected objects. The 
small aperture was necessary because of crowding. 

The near-infrared data used in this paper were de- 
rived from the instrumental magnitudes obtained in the 
non-standard filter systems of NICMOS, using the infor- 
mation of the mean metallicity based on optical colors. 
The relation between (V — I) and the metallicity is rather 
well calibrated, using the globular clusters of the Milky 
Way or M31, respectively. However, considering the fact 
that the correction between F160W and the Cousin EE- 
band for a mean metallicity of (Fe/H) —0.3 (IC 4051) 
and (Fe/H) —0.75 (NGC 3311), as described in Sect. l2~2~Tl 
and 12. 2. 31 differ only insignificantly (0.115 mag in IC 4051 
compared to 0.110 mag in NGC 3311) leads us to the con- 
clusion that the a priori assumption of the metallicity does 
not alter the final outcome of our analysis. Nevertheless, 
since additional near-infrared data are desirable for an im- 
proved GC sample in both galaxies, observations using 
broad band filters with a response function closer to the 
Johnson-Cousin system are clearly preferable, or alterna- 
tively observations in multiple bands, allowing a proper 
color correction. 

2.2. Near-Infrared Data: HST/NICMOS2 

Both galaxies were observed with NICMOS Camera 2 
(fov: 1972x1972, 07075 pix" 1 ) through the F160W filter 
(approximately H band) in a slightly overlapping (^10-20 



pixels) four-point mosaic pattern, centered on the galaxy. 
The combined mosaic image covers approximately the 
same area and galaxy field as the WFPC2 PCI chip. Each 
of the four quadrants is a combination of 5 sub-images. 
The exposure time for each sub-image is 512 sec, the 
total exposure time sums up to 2560 sec (see Table EJl- As 
for the optical filter bands, we give the details for data 
reduction separately for each target. 



2.2.1. IC 4051 

In IC 4051 the SExtractor routine 
HBertin fc Arnout 199 6:) was used to identify all tar- 
gets in the four quadrant images of the galaxy-subtracted 
NICMOS2 image. This resulted in the detection of a 
total of 642 targets. The aperture photometry was then 
carried out using the non-standard IRAF task CCDCAP 
(Mig hell fc Rich 19 95), which is optimized for aperture 
photometry of undersampled images. Hereafter, the pho- 
tometry procedure recommended in the NICMOS Data 
Handbook 1 (Dickin son et al. 2002j) was applied. The 
instrumental magnitudes were measured with an aperture 
of 075 (radius, corresponding to 6.6 pixel) and corrected 
to "infinite" aperture radius by adding -0.1517 mag, 
calibrated by adding 23.566 mag (which is the magnitude 
of a star with count rate of 1 electron per second in 
the NICMOS F160W filter), and corrected for the total 
quadrant image exposure time. Finally, we performed a 
transformation from the NICMOS F160W filter magni- 
tudes to standard H magnitudes ( |Stephens et a l. 2000). 
Based on the mean metallicity of IC 4051 globular 
clusters of (Fe/H) —0.1 HWoodworth & Harris 2000|) . 
combined with observations of Galactic and M31 globular 
clusters ( |Barmby et al. 2000) ) , leads to the following 
transformation: 

H = F160W — 0.115 (±0.048) mag (2) 

After all calibration corrections were applied to the 
instrumental magnitudes, a number of targets were 
rejected: all targets within 12 pixels of any edge of a 
quadrant image, all targets for which CCDCAP failed 
to measure a magnitude, all targets with SExtractor 
star/galaxy indices less than 0.5 (indicating they are 
likely galaxies), targets located near the the center of 
the subtracted galaxy (within ~270, where photometric 
errors are very large due to the background galaxy light), 
and for identical targets located within the overlapping 
regions of two quadrant images, the one with the larger 
uncertainty. A total of 348 objects remained. After trans- 
forming the WFPC2 pixel coordinates to the NICMOS 
coordinates, and allowing an offset of up to 2 pixels in 
both x and y coordinates, the optical and near-infrared 
data were matched. The final sample now includes 256 
objects. 



1 See http:www.stsci.edu/hst/nicmos/ 
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2.2.2. NGC 3311 

As in the optical filter bands the source detection on the 
NICMOS images was done using DAOPHOT/PHOT as 
well as SExtractor. To be as consistent in the data reduc- 
tion as possible, we will stay with the SExtractor results. 
The transformation of the FI60W data to the Johnson 
H band followed the procedure described in Section 12. 2. II 
and results in the following correction: 

H = F160W - 0.110 (±0.024) mag (3) 

The differences in the transformation between 
F160W and H-band magnitudes between NGC 3311 
and IC 4051 originates in the difference of the mean 
metallicity, derived from the (V — I) distribution, where 
we have adopted a mean metallicity of ( Fe/H )~-0.75 
(Bro die, Larsen, fc Kissler-Patig 2000| ). Besides the 
F160W- H band transformation, the instrumental magni- 
tudes were aperture corrected, the photometric zero point 
for NICMOS2 applied and finally corrected for galactic 
reddening (see Table After combining optical and 
near-infrared data the globular cluster sample includes 
148 objects. 

In this series of papers we concentrate on the age 
structure in globular cluster systems, which we derive 
from color-color diagrams. The basic tool in this approach 
(see Section l3~2l and |Hempel fc Kissler-Patig 20 04 ) is the 
color distribution of the observed globular clusters with re- 
spect to SSP model isochrones. Therefore we have to con- 
sider the contamination of the globular cluster samples. 
Considering the color range, the contaminants are most 
likely background galaxies. Although we apply selection 
criteria, e.g. shape parameters and color limits, we are 
aware that contamination is still an issue. Nevertheless, 
the final globular cluster samples, combining optical and 
near-infrared data, contain only objects in the very cen- 
tral region of the CCS, where we find the highest globular 
cluster density. In combination with the sample size («100 
objects), as discussed in Section l3~2l we estimate that un- 
resolved background galaxies make only a minor contri- 
bution to our sample and should have a minimal effect on 
the age distribution. 

3. Results 

3.1. Color-Color Distributions 

Our basic tools to derive ages and metallicities of 
stellar populations from photometric data are color- 
color diagrams together with various SSP models; for 
example by Bruzual & Chariot (2000, 2003), Vazdekis 
(1999), Maraston (2001) or their most recent release 
UMaraston 2005|) . We will use the Bruzual & Chariot 
models IjBruzual fc Chariot 2003| throughout the main 
body of the paper whereas Appendix ^] is based on the 
Maraston models IjMaraston 2005J) which indicates that 
only minor differences occur between models for ages >1 



Gyr and that our main goal of detecting different aged 
subgroups is not compromised. In our further analysis 
we select only clusters which satisfy an error limit of 
fi(V-i) an d S(v-H) < 0.15 mag. The error limit was 
derived as a combination of the photometric error in 
two filter bands with <5<0.1mag each, corresponding to a 
S/N=10 for detected objects. The error selected globular 
cluster samples for NGC 3311 and IC 4051 contain 105 
and 149 objects, respectively. In our further analysis of 
the age structure additional color cuts will be applied 
to minimize contamination with unresolved background 
objects and to select objects within a restricted color 
range that allows the best age resolution (see Section 
13.2(1 . With respect to previous studies, e.g. by Hilker 
(NGC 3311, Hil ker 20021 IHilker 2003a) and W oodworth 
& Harris (IC 4051, IWoodworth fc Harris 2000)> we note 
the much smaller size of our sample, which is mostly due 
to the included H-band data, and the small field of view 
of the NICMOS2 instrument. As we will discuss later 
(see Sect. 13.21 and IBT). the limited spatial distribution of 
the globular cluster sample becomes of importance for 
our interpretation of the derived relative sizes of the age 
sub-populations. 



3.1.1. IC 4051 

In Figures [3 and |31 we show the color-magnitude diagram 
(CMD) and the (V — I) vs. (V — H) color-color diagram 
for IC 4051. Clearly, the data below the 50% completness 
limit of the H photometry are unreliable and will be 
excluded from further analysis, as will all objects with 
color errors above our error cut. In Woodworth & Harris 
(2000) it was tentatively suggested that the (V — I) 
distribution of the IC 4051 globular clusters is bimodal 
and, therefore, we attempted to fit a double Gaussian 
function to the {V — if )-color distribution, which is three 
times more sensitive to metallicity than (V — I). If the 
complete sample set is included in this analysis, then 
comparably good fits are obtained for variable amplitudes 
for both sub-populations or (forced) equal amplitude, 
equal width settings. In the latter case the Gaussian 
centers are at (V — H)=2.70 and 3.18. The mean of these 
two Gaussian centers is 2.94, which is approximately 
equal to both the mean of all colors and the center of the 
single Gaussian fit. Thus, this fit essentially devolves to 
merely reproducing the single Gaussian case. However, 
since uni- and bimodal fits give similarly good results 
(see Figure 0J, we will stay with the assumption of a 
unimodal (V — H) color distribution and a corresponding 
mean color of (V — H) = 2.92. We note that the situation 
changes drastically if we apply the selection by error 
limit. As we can see in Figure |21 the blue population 
(V — H <2.4) is significantly diminished by the error 
limit, which makes it difficult to draw strong conclusions 
from the color distribution. The metallicity distribution 



6 



Hempel et al.: Extragalactic Globular Clusters in the Near Infrared V 




1 2 3 4 5 

V H |mag] 



Fig. 2. V vs. (V — H) color-magnitude diagram for IC 4051. The top sub-panel shows the color distribution of all (open 
histogram) and error-selected (hashed histogram) objects. The lower sub-panel shows the color-magnitude diagram. 
Herein the filled symbols mark the selected clusters, while the open circles mark rejected objects. The dashed line 
marks the 50% completeness limit for the H-band. 

T~> 1 1 1 I ' 1 1 ' I 1 1 1 1 I 1 1 1 1 I 1 ' 1 1 

1£ : IC4051 t ; 




5 Gyr / 1 Gyr 



0.6 - 1 1 



1.5 2 2.5 3 3.5 4 

V-H [mag] 

Fig. 3. (V— I) vs. (V—H) color-color diagram for IC 4051. All data are corrected for foreground reddening (see Table^). 
The 15 Gyr, 5 Gyr and 1 Gyr isochrones are based on the Bruzual & Chariot SSP models i|Bruzual & Chariot 2 003). 
The solid squares mark the metallicity in [Fe/H], increasing with (V — H) from -2.25 up to +0.56 . The box marks 
the color range selected for the determination of the cumulative age distribution. 



is further investigated in Section 

With the discovery of bimodal color dis- 
tributions in globular cluster systems (e.g. 
IWhitmore fc Schweizer T995I Geisler, Lee & Kim 1996; 
Kundu & Whitmore 2001a; Gebhardt & Kissler-Patig 
1999) the interest in the GCSs of early-type galaxies, 
which were supposedly old and uniform, grew immensely. 



The spatial distribution of GCs adds an extra clue to 
their origin. For example, in the Galaxy, metal rich GCs 
are found to be more confined to the Galactic disk out 
to distances up to ~10 kpc from the Galactic center 
( Sca rle & Zinn 1978j) . and metal-poor clusters are found 
at even larger distances but spherically distributed in 
the Galactic halo (Ash man fc Zepf 1998D . Although the 
near-infrared data for IC 4051 cover only the innermost 



Hempel et al.: Extragalactic Globular Clusters in the Near Infrared V 
Table 3. Mean Color versus Radius for IC 4051 
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30 


2.901 


0.745 


5.0-7.5 


28 


2.897 


0.530 


7.5-10.0 


41 


2.938 


0.478 


10.0-12.5 


51 


2.737 


0.635 


12.5-15.0 


51 


2.984 


0.612 


15.0-17.5 


38 


3.044 


0.724 


17.5-20.0 


25 


3.114 


0.451 



region of IC 4051 (maximum distance ^9 kpc) it is still 
worth investigating the radial distribution of the globular 
clusters mean colors. Figure shows the distribution of 
globular cluster colors as a function of galactocentric 
angular distance R; we have listed the mean colors in 
2'.'5 wide bins in Table |3| We note that there is no 
obvious trend in the mean colors. V — H is primarily 
sensitive to metallicity, so this indicates the lack of any 
strong metallicity gradient in the IC 4051 GCS, for radial 
distances R<22"5. 




(V - H) 

Fig. 4. (V — H) color distribution in the complete IC 4051 
sample showing a single Gaussian fit (middle) and two dif- 
ferent double Gaussian fits (top: equal amplitude/equal 
width, bottom: variable amplitude/ variable width). The 
dotted lines mark the contribution of each individual com- 
ponent, solid lines are their sum. 

Comparing the combined optical and near-infrared col- 
ors with SSP model isochrones (see Figure [3J reveals that 
only a fraction of the red cluster population is consistent 
with an old (>10 Gyr) population. We find a large frac- 
tion of objects whose (V — I) color at a given (V — H) 
is too blue to fit into an "old" globular cluster scenario. 




0.5 1 

distance fr_effl 



1.5 



Fig. 5. Distribution of globular cluster colors for IC 4051 
as a function of radius from the center of the galaxy (open 
circles) . The dashed line marks the mean value of the en- 
tire data set, (V-H}=2.922 mag. The solid circles show the 
mean colors in increasing radius bins, as listed in Table [3J 

Their color is much better matched by an intermediate 
age population (i.e. with an age between ~1 and 7 Gyr). 
Such a large intermediate/young population of globular 
clusters seems suspicious, since the integrated galaxy light 
does not show any hint of a recent major star formation 
event and it is a relatively isolated galaxy. With respect 
to the integrated light argument we refer to the work by 
Puzia et al. (2002) on NGC 4365 who found a similar sce- 
nario where the combined (V — I) and (V — K) color of 
the host galaxy is consistent with a metal-rich and rela- 
tively old stellar population although the globular cluster 
system indicates a significant fraction of intermediate age 
globular clusters as well as old ones IjPuzia et al. 2 002 
|Hempel fc Kissler-Patig 2004| ). The apparent contradic- 
tion between the two results breaks down if we compare 
the observed fraction of globular clusters with the stel- 
lar content of the galaxy. The derived ratios between age 
sub-populations in both galaxies refers to only a small 
sub-sample of the globular cluster system and is therefore 
subject to various bias effects, as will be discussed later. 
In this respect we also note the most recently published 
results on NGC 4365 by Brodie et al (2005). In this paper 
spectral data of GCs are used to derive their ages, which 
seem to contradict previously obtained age estimates, both 
spectroscopic and photometric, and claim an average age 
of 11 Gyr for the complete GC sample, with no evidence 
for intermediate age clusters. At this point we have to take 
into account that the different samples have only few glob- 
ular clusters in common. However, the authors discuss the 
possibility of spatial distribution effects, which is indeed 
an important issue when drawing conlusions on ages of se- 
lected samples. Clearly more data are required to explain 
the different results. 

At this point we want to emphasize that only within 
the red sub-population is it possible to distinguish the 
age of globular clusters whereas the blue population 
should represent only the first generation, but deriving 
even crude age estimates for such blue clusters is very 
difficult. This supposedly old fraction is NOT included 
in the derived size of the different age populations. In 
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Fig. 6. V vs.(V — H) color-magnitude diagram for NGC 3311. The top sub-panel shows the color distribution of 
all (open histogram) and error-selected (hashed histogram) objects. The lower sub-panel shows the color-magnitude 
diagram. Herein the filled symbols mark the selected clusters, while the open circles mark rejected objects (by error 
cut). The dashed line marks the 50% completeness limit for the H-band, derived from the detection rate in artificial 
data sets using the Addstar routine in the IRAF package. 
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Fig. 7. (V — I) vs. (V — H) color-color diagram for NGC 3311. The data obey the selection by photometric error 
cut (<0.15 mag) and are corrected for foreground reddening (see Tabled]). As examples the 15 Gyr, 5 Gyr, and 1 
Gyr isochrones, corresponding to the Bruzual & Chariot model isochroncs (IBruzual fc Chariot 20 03) are marked. The 
metallicity, rising from [Fe/H]=-2.25 to [Fc/H]=+0.56, is marked by solid squares. The color range selected for the 
determination of the cumulative age distribution is marked by the box. Compared to IC 4051 the (V — H) color range 
is slightly shorter to fit the observed data. 



Section ^ we mentioned a gas rich merger as a possible 
scenario to form a second generation of globular clusters, 
which requires a certain minimum galaxy density. IC 4051 
is found somewhat outside the core of the Coma galaxy 
cluster, surrounded by numerous other galaxies, but not 



in the center of the gravitational potential of the galaxy 
cluster. However, additional evidence for a less passive 
evolution of this galaxy is the existence of a co-rotating 
galaxy core, which rotates at a different velocity from 
the main galaxy body. It is not known yet whether it is 
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dynamically decoupled or not (M ehlert et al. 1998|l . Such 
peculiarities are quite common in early-type galaxies (e.g. 
IDavies et al. 20011 Ide Zeeuw et al. 2002(1 and we do not 
intend to link the existence of such a co-rotating core 
directly with a second star formation event. However, if 
the dynamical structure of an early-type galaxy sets it 
apart from the formerly- uniform and simple structured 
picture, then these will be prime targets to look for 
multimodal age distributions. In this respect we want 
to emphasize that our classification of IC 4051 as an 
early-type galaxy in a less dense environment is rather 
relative, since the Coma galaxy cluster represents a high 
density environment (projected central galaxy density 
of «39 galaxies/Mpc 2 , IKent fc Gunn 1982|l compared 
to the Hydra I galaxy cluster (17 galaxies/Mpc~2, 
Richtcr 1989). Nevertheless, due to its distance of 15f5 
(~460 kpc, see |Jorgensen fc Franx 199"4| ) from the cluster 
center we do not consider IC 4051 as a galaxy in a high 
density environment anymore. 



3.1.2. NGC 3311 

The V vs. (V — H) color-magnitude diagram and (V — 
I) vs. (V — H) color-color diagram for NGC 3311 are pre- 
sented in Figures|H|andEI The histogram in the upper part 
of the CMD shows the color distribution of the complete 
sample as an open diagram, whereas the error selected 
sample is shown by the hatched histogram. 

The color-color distribution in Figure [7| shows that 
in addition to many old GCs, a large fraction of red 
globular clusters (V — H > 2.5) have a (V — I) color 
which is too blue to be consistent with an old cluster 
population. A direct comparison to the SSP isochrones 
assigns those objects an age <5 Gyr. The possibility 
of a young/intermediate population of globular clusters 
formed in a merger was already discussed by Hilker 
(2002). Our near- infrared observations do not favor the 
second possibility, viz. an old and metal poor population. 
Interestingly, NGC 3311 contains a prominent dust lane 
in its center (e.g. |Vasterberg, Jorsater fc Lindblad 1991| 
Grillmai r et al. 1994|l . which adds further evidence to 
support a recent merging or accretion event. The dust 
content in the core of NGC 3311 was estimated as 
10 3 - 5 M Q IjFerrari et al. 1999jl based on V and R-band 
photometry. Grillmair et al. (1994) derived a lower limit of 
4.6xl0 4 Mq for the total dust mass. Assuming a similar 
gas-to-dust ratio as found in the Galaxy, wigcs/indust™ 

100 

IjJura 19861 |McNamara, Bregman fc O'Connell 1990| ), 
the total gas and dust mass will reach at least 4.6 x lO e M0 . 
The globular cluster sample included in our analysis, how- 
ever, refers only to the innermost region in NGC 3311. 
Additional near-infrared observations covering larger 
galactocentric radius will be necessary before any con- 
clusion can be inferred from comparing our small central 
sample to the entire blue globular cluster population 
found by Hilker (2002), using optical photometry alone. 



As previous studies on NGC 4365 l|Puzia et al. 2 002 ) 
have shown, results based on optical photometry alone do 
not always indicate the existence of age sub-populations. 
The distribution of (V — I) colors shows only the common 
bi-modal distribution, mostly caused by the metallicity. 
Only the combination with near-infrared colors reveals the 
more complex age structure of the globular cluster system. 



3.2. Cumulative Age Distribution 

Our technique to derive cumulative age distributions, 
as described in Hempel et al. (2003) and Hempel & 
Kissler-Patig (2004), is to help lift the age- metallicity 
degeneracy ( |Worthey 1994| ) via the combination of optical 
and near-infrared broad band colors. Hereby the higher 
sensitivity of infrared colors to metallicity as compared 
to age is utilized. An example of the comparison of the 
age and metallicity dependency of (V — I) and (V — K) 
(~ (V - H)) can be found in Puzia et al. (2002) and 
references therein. The cumulative age distribution is 
then defined as the percentage of objects older than T 
Gyr, and is derived by comparing their (V — I) color 
with the SSP model prediction for a population that 
is T Gyr old. The lowest age bin is assigned to Gyr 
and includes all GCs within the selected (V — H) color 
range. With respect to a possible formation scenario of 
the GCSs and their host galaxies, we want to go one step 
further and set constraints on the age and the relative 
size of the cluster sub-populations. The method is based 
on the comparison between observed and simulated 
GCSs Pempel et al. 2003| |HenTpel fc Kissler-Patig 2004} 
via a x A -test. To simulate the color distribution in 
GCSs and to derive the age distribution in simulated 
and observed systems, the SSP model isochrones of 
Bruzual & Chariot (2003) were again used. Examples 
for the simulated color-color diagrams and the resulting 
cumulative age distributions are shown in Figure [HJ 
Instead of a step-function-like age distribution we find a 
much more gradual slope in the age distribution, which 
is mostly due to the fact that the (V — /) colors in our 
simulations are not strictly independent from metallic- 
ity effects and, therefore, show a certain dependency 
on (V — H). The final consequence is that even with 
a much improved age resolution of combined optical 
and near-infrared colors, the determination of accurate 
absolute ages for single globular clusters is still out of 
reach. Our intention is to detect age sub-populations, to 
set constraints on their age and to find systematic effects 
of external parameters on the age structure. In nearby 
galaxies more precise age determinations can be derived 
from spectroscopic spot checks. For our simulations we 
assume a composition of two age sub-populations with 
ages ranging from 1, 1.5, 2, 3, 5, 7, 10 and 13 Gyr. 
Hereby we vary the amount of young/intermediate age 
globular clusters between 100% and 0%. The x 2 ~ test 
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Fig. 8. Simulated color-color diagrams (left) and the resulting cumulative age distribution (right) for various age 
combinations. The two populations are equally mixed between the old and intermediate age populations. The ideal 
colors, following the Bruzual & Chariot (2003) SSP models, were smeared randomly with an up to 3tr photometric 
error. 



finds the best fitting model with a given age and size ratio. 

As shown in Sect. 13. Ol and 15.1. li the color-color dia- 
grams give us a first indication for the existence of multiple 
age sub-populations in the GCSs of both NGC 3311 and 
IC 4051. If we compare the cumulative age distributions 
(Figure El with the simulation of a purely old system (see 
Figure |H1 upper, right panel) the differences are obvious. 

From Figure the age distribution in both GC 
samples assigns ^65% of the globular clusters to an age 
younger than 5 Gyr. Recall that this refers only to GCs 
with V — H > 2.5. Although the combination of optical 
and near-infrared colors helps lift the age-metallicity de- 
generacy, (V — I) is still mildly affected by the metallicity. 
Combined with the scatter in both colors, induced by the 
photometric errors, and the relatively small sample size 
(IC 4051: 101 GCs, NGC 3311: 70 GCs), this results in a 
significant uncertainty with respect to the size estimate 
of the sub-populations. To overcome this obstacle we 
will not use the cumulative age distribution as a direct 
measure for the age structure in the GCSs, but rather 
compare the observations with the simulated systems 
via a x 2 - test. In this test we determine the difference 
between the cumulative age distribution in observed 
and simulated systems. Nevertheless, we still have to 
accept a certain degeneracy between age and size of 
sub-populations, which results in a finite accuracy with 
respect to the best fitting model. Since the major goal 
of this project is to detect age sub-populations and set 
some constraints on their age, it is of no importance here 
(and impossible for us to detect) whether the second 
generation of globular clusters is 2 or 3 Gyr old and 
contains 30 or 40 % of the globular cluster sample. 
This uncertainty becomes even more important for 
NGC 3311, where a large fraction of the observed GCs 
are excluded from the age analysis by the color cut. In 
the comparison of our results with previously obtained 
age estimates we will therefore concentrate on IC 4051 



(e.g. IMehlert et al. 19981 IMehlert et al. 2003|l. Base d on 
spectroscopic data (Figure 8 in IMehlert et al. 1 998 ) the 
age of the galaxy core of IC 4051 was determined as 12-17 
Gyr, based on the Worthey SSP models jWorthey 1994| ). 
We note that agreement with model predictions was 
only obtained for the innermost 5'.'0 of the galaxy core, 
for which the photometric data of the globular clusters 
suffer from contamination by the diffuse galaxy light. 
An additional error source is the fact that in the diffuse 
galaxy light the contributions of stellar sub-populations 
are blended and hard to disentangle. 

The results of these x 2 tests are presented in Figures 
ITU1 and ^2 The \ 2 test results find the best fitting 
model to the age distribution in both GCSs to be a 
composite of an old population and a second popu- 
lation which is ^10 Gyr younger. In both galaxies, 
the best fitting model to the age distribution suggests 
approximately 60% to 70% intermediate age clusters. 
Such a surprisingly large fraction of young/intermediate 
globular clusters, also found in NGC 4365 and NGC 5846 
(He mpel fc Kissler-Patig 2004| , needs confirmation and 
careful investigation, especially since such large fractions 
of young stars should be visible in the integrated light, but 
are not detected yet. In fact, previous investigations based 
on optical photometry IjWoodworth fc Harris 2000jl . sug- 
gest the bulk of the GCS being formed in situ, which, in 
the authors understanding, also includes the possibility 
of merging of almost completely gaseous objects. A late 
gaseous merger is not in any way ruled out using as well 
the argument of the central co-rotating disk and its very 
high metallicity. 

The major goal in this work is to apply a semi- 
numerical method for detecting globular cluster sub- 
populations based on optical and near-infrared colors. As 
described in Hempel et al. (2003) and Hempel & Kissler- 
Patig (2004) this procedure has been developed in order to 
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Fig. 9. Cumulative age distribution in NGC 3311 (open squares) and IC 4051 (solid circles) derived following the 
procedure given in Section r3~2l The left panel shows the absolute numbers of globular clusters older than a given age. 
In the right panel the distribution was normalised with respect to the total number of objects found. If we assume 
the blue population of globular clusters ((V — H)< 2.5) to consist of only 13 Gyr old objects the cumulative age 
distribution would change as shown by the grey patterns. 



partly compensate for the photometric uncertainties, in- 
creased by the usage of near-infrared observations, when 
deriving age constraints. Instead of a direct comparison 
between the colors of each individual observed globular 
cluster with model predictions, as attempted with color- 
color diagrams (see Figures |3] and we rather look at the 
density of GCs with respect to a set of model isochrones 
and compare it to simulated GCSs. Although the results 
of both approaches may differ in the actual age estimates, 
the existence of distinct age sub-populations should be re- 
vealed by both approaches, the semi-numerical one being 
less sensitive to photometric errors. Again, concentrating 
on the IC 4051 results shown in the color-color diagram 
(see Figure 2J) we are unable to derive individual globu- 
lar cluster ages by comparing the GC colors directly with 
the SSP model predictions, but nevertheless find a large 
fraction of globular clusters with (V — I) colors being bet- 
ter represented by a isochrone for a 5 Gyr population or 
even younger. A better agreement between the age esti- 
mate via our semi-numerical approach should be possible 
by refining the age grid for the Monte-Carlo simulations 
(e.g. |Hempel fc Kissler-Patig 2004| | as well as for the de- 
termination of the cumulative age distribution. 

The derived age structure will therefore be affected by 
different factors, the most important being stated below: 

— Firstly, we have to consider that the sample used 
in our analysis is color-selected, so only objects with 
2.5 < (V - H) < 3.5 and 0.6 < (V — I) < 1.4 are 
used. The lower limit to (V — H) is set due to the 
small color differences in the SSP models. The blue, 
metal-poor population of globular clusters will mostly 
contribute to the old sub-population and change the 
size ratio between old and young population in favor 
of the first generation objects. The upper limit, as 
well as the (V — I) interval, is set to avoid signifi- 
cant contamination of our sample with background 
galaxies (P uzia et al. 2002|) . Especially the (V — H) 
blue fraction, supposedly old globular clusters, will 
affect this result and shift the ratio between old and 
young/intermediate age globular clusters to higher 



values. However, as shown in Figure even in the 
extreme case where all GCs bluer than (V — 7J)=2.5 
are old, 13 Gyr old clusters there remains a substantial 
fraction of young/intermediate age clusters. 

— Secondly, we have to be aware of spatial bias effects. 
Due to the small NICMOS field of view, centered 
on the galaxy, the data cover only the innermost 
region of NGC 3311 and IC 4051. As observed by 
Schweizer & Seitzer (1998) and simulated by Hibbard 
& Mihos (1995) for the merger remnant NGC 7252, 
the central region of a galaxy is where one would 
expect to preferentially find the second generation of 
globular clusters, if they are indeed formed during 
a merger or accretion event. Although mostly born 
in the tidal tails ( |Li, Mac Low fc Klessen 2004| ) of 
a merger, these star clusters will funnel toward the 
center of the merger remnant within several 100 Myr 
up to 1 Gyr (H ibbard fc Mihos 1995(1 . Our samples 
are therefore spatially biased and the fraction of 
young/intermediate age globulars should be overesti- 
mated, compared to the global fraction over the entire 
galaxy. 

— Thirdly, small number statistics will affect the % 2 - 
test. As shown in Hempel & Kissler-Patig (2004) the 
results of this test become unstable (i.e. show multi- 
ple best fitting age/size combinations) for small data 
sets. Previous tests have shown that >100 objects are 
needed to obtain reliable results. Strictly speaking this 
condition is only (barely) achieved in the IC 4051 
sample. Recall that the total sample of IC 4051 and 
NGC 3311 contain 256 and 148 objects, respectively, 
but after applying the selection criteria (photomet- 
ric error and color cuts) only 101 (IC 4051) and 70 
(NGC 3311) objects contribute to the age distribu- 
tions. Although the presence of a significant popu- 
lation of young/intermediate age globular clusters in 
NGC 3311 agrees with various other features of the 
galaxy (e.g. central dust lane) more data are required 
for a solid argument. 
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Fig. 10. x 2 - test result for NGC 3311. Different levels represent the reduced x 2 of the comparison between the 
cumulative age distribution in NGC 3311 and various models. The best fitting model consists of ~ 60% intermediate 
age globular clusters and ^40% of 13 Gyr old objects. The large size of this young population is surprising and is 
discussed in Sectional 
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Fig. 11. The comparison between the age distribution in IC 4051 and the simulated systems again (as in Figure HU|l 
reveals a large fraction of young/intermediate age objects (>~60%), with the population of second generation globular 
clusters slightly larger than in NGC 3311. 



— Fourthly, as mentioned in Section ^ the SSP model 
isochrones by Bruzual & Chariot (2003) have been ap- 
plied. As shown by Hempel et al.(2003) and Hempel & 
Kissler-Patig (2004) the color predictions for a given 
age by various SSP models differ particularly in the 
red color range and model dependencies should be 
considered and wc will do so in detail in Appendix 
lAl Since our method is based on the object density 
in color-color diagrams with respect to a given set of 
SSP isochrones, the final result, i.e. the detection of 
age sub-populations will hold independently of the SSP 
model. The number of objects required for a stable re- 
sult of the x 2- test is however strognly dependent on 
the 8a.ge/5(V — I) gradients. The models by Maraston 
(2005) are somewhat hampered by the fact that stellar 
populations with ages below 5 Gyr (below 2 Gyr for the 
2005 release) show a non-monotonic increase of (V — I) 
with (V — H) relation, which leads to confusion in our 
counting algorithm. The major advance in the more re- 
cent SSP models are the inclusion of a well developed 
AGB phase (e.g. Renzini 1981, Maraston et al.2001b), 
the possibility of a different [a/Fe] -enhancement 



({Thomas &: Maraston 2003|l and the structure of the 
horizontal branch (Maraston 2005), resulting in color 
shifts. Although a- enhancement seems to be more im- 
portant when objects with solar metallicity are inves- 
tigated no solid conclusions can be drawn, due to less 
solid calibration of the models. We will include the 
latest release of the Maraston models (2005), target- 
ing the influence of the horizontal branch structure in 
an additional section (see Appendix E} ■ 



4. Metallicity Distribution 

In the framework of galaxy formation and evolution stud- 
ies the age structure of the galaxy and its GCS plays a ma- 
jor role hence our main interest is focused on the age struc- 
ture. Nevertheless, since the latter is determined by using 
SSP model isochrones, which combine age and metallicity 
predictions, it is necessary to test to what degree age and 
metallicity effects are still entangled. In the following we 
present the results of rough metallicity estimates based on 
the derived age structure. 
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A 2 Gyr old stellar population, with a (V—H) color be- 
tween 2.5 < (V — H) < 3.5, i.e. consistent with the color 
limits set for the determination of the age distribution, 
contains objects with metallicities between [Fe/H]=+0.00 
and [Fe/H]=+l.l, respectively, as derived via linear inter- 
polation using the SSP models. Assuming a mean color 
of (V - iJ)=3.0 we obtain [Fe/H]=+0.56. Although, as 
shown in Figure ^2 the best fitting model consists of a 
2 Gyr old second globular cluster population, the age un- 
certainty allows as well a higher age, i.e. ~5 Gyr. In this 
case the metallicity of objects in the same color range 
lies between [Fe/H]=-0.2 and +0.9, with a mean value of 
[Fe/H] =+0.4 for (V-H)=3.Q. As we will show, this corre- 
sponds closely to the mean metallicities determined for the 
young/intermediate GCs in IC 4051. Since the metallic- 
ity depends quite strongly on the derived age (more than 
the other way around) and although we can detect age 
sub-populations and set constraints on the age, the age 
uncertainty allows a large variety of metallicities. The fol- 
lowing discussion about the metallicity distribution has to 
be seen in the light of this persistent degeneracy between 
age and metallicity. 

With the introduction of combined optical and near- 
infrared observations we are able to investigate the metal- 
licity distribution approximately independently from the 
age structure, with the above caveats. In order to derive 
globular cluster metallicities we use the above-derived in- 
formation about the age structure in the observed glob- 
ular cluster sample (see Figures and . The results 
in Section l3~2l tell us that the age structure in both sys- 
tems is well matched by a mix of a 13 Gyr old population 
with globular clusters which are ~10 Gyr younger. The 
determination of the metallicity distribution follows the 
procedure given below and applies the Bruzual & Chariot 
(2003) SSP model isochrones. 

— We split the globular cluster sample (selected by pho- 
tometric error) with respect to the (V — H) color into 
a blue ((V — H)< 2.5) and presumably old population 
and a second, red (2.5 < (V — H)) sub-population, con- 
taining both old as well as young/intermediate glob- 
ular clusters. In addition we set an upper limit of 
(V — H)< 3.5 to avoid strong contamination with back- 
ground galaxies ( |Hempel fc Kissler-Patig 2004| ). 

— For globular clusters found in the blue sub-population, 
the metallicities (given as Z/Z©) were calculated by 
interpolating the metallicity- (V — H) correlation of the 
13 Gyr model isochronc. 

— Corresponding to the age structure, we found that 
the red sub-population is split into old and interme- 
diate/young populations, and use the 5 Gyr model 
isochrone for the latter. All objects with (V—I) < (V — 
1)5 Gyr are then assumed to form the second genera- 
tion of globular clusters. Although the comparison of 
the cumulative age distributions to simulated GCSs as- 
signs the intermediate/young objects to an age of ~2 
Gyr, we allow an age uncertainty of 3 Gyr, based on 
the x 2 contours, and set the age split at 5 Gyr. 



— The metallicity is then determined by interpolating 
Z/Z© at (V - H) GC using the 13 Gyr and 2 Gyr 
isochrones (Figures 1121 and IT^l left panels). For com- 
parison we also derive the metallicity distributions 
assuming a purely old globular cluster population 
(FigureslT^land lTSI middle panels) and for mixed GCSs 
containing 13 Gyr and 5 Gyr old objects (Figures ITS! 
and FLU right panels). In both Figures the upper pan- 
els give the metallicities on a linear scale, whereas the 
lower panels refer to [Fe/H] values. 

— The metallicity estimates depend strongly on the age 
structure and are therefore also hampered by bias ef- 
fects (see Section l3~2"|) . The case of NGC 3311 suffers 
from small-number statistics, which makes the results 
of the x 2- test statistically unstable. Nevertheless, the 
comparison of the (V — H) colors, seen in Figures El 
and [7| shows that NGC 3311 lacks the most metal rich 
globular clusters found in IC 4051. Due to the above 
uncertainty for NGC 3311 we will concentrate on the 
GCS in IC 4051. 

In the discussion about the bimodality of metallicity 
distributions we want to emphasize the importance of the 
chosen scaling. It is interesting to note that the metal- 
licity distribution for a purely old (13 Gyr) population 
seems to suggest a bimodality in the [Fe/H] distribution 
for the IC 4051 GC, which is NOT seen if we transform 
the GC metallicities into linear abundances. However, if 
we fit the [Fe/H] distribution by either a double or a sin- 
gle Gaussian, as shown in the middle panels in Figure IT4l 
the resulting \ 2 of the fit is found to be x 2 =0.6 and 0.8 
for the double and single Gaussian, respectively, virtually 
identical. Hence there is no strong evidence for bimodality 
in the metallicity distribution, which considers all clusters 
as old. Although a bimodal distribution is not favored over 
unimodality in the 13 Gyr old population, it becomes ob- 
vious in the samples assuming mixed globular cluster age 
distributions (see upper and lower panel in Figure [14^1 . 

Since wc find objects within a wide range of metallic- 
ities, applying a logarithmic scaling, as given in the dis- 
tribution of the [Fe/H] values, increases the chances for 
detecting metallicity populations. In the linear Z space 
formerly clearly detected [Fe/H]-populations are much 
more spread out and may not be detected anymore, 
e.g. in smaller samples or in globular cluster populations 
with a smaller age difference. The metallicity distribu- 
tion given in Z/Z Q for a 13 + 2 or 5 Gyr old popula- 
tion is widely spread. We caution that a bimodal color 
distribution does not necessarily translate into a clear bi- 
modal metallicity or linear abundance distribution (see 
also |Kissler-Patig 2002D . 

The mean metallicity based on the (V — I) 
colors alone and a unimodal color distribution 
QWoodwortfa fc Harris 2000) was determined to be 
(Fe/H)=-0.3 or Z/Z =O.5. Working with the more metal 
sensitive (V — H) color we obtain the following results. 
Assuming a purely 13 Gyr old population in IC 4051 (see 
Figure 1121 central panel), we derive a mean metallicity 
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Fig. 12. Metallicity distribution in IC 4051 assuming different GC distributions. We show the metallicity distribution 
both in terms of linear metal abundance (upper panel) as well as in logarithmic metallicity [Fe/H] (lower panels). In 
the left panels the metallicity distribution is based on the assumption of a 2 and 13 Gyr old population (following the 
best fit model from the age distribution), whereas the central panel assumes all clusters to be 13 Gyr old. We also 
show the metallicity distribution assuming the young/intermediate globular clusters to be 5 Gyr old (right panel). The 
color coding of the histograms is as follows; dark grey: young and metal rich, light grey: old and metal rich, black: old 
and metal-poor. 



40 

u 30 

u 

g 20 

eJ 

Pi 



10 


30 



" I " 1 1 1 1 1 " 1 1 
13+2 Gyr 



■fisKiNFbi iYi 1j i\h= 



0.5 1 1.5 2 
z/ solar 



v 20 

a 



I " " I " " I '. 

13+2 Gyr: 




-2-10 1 
z/ solar 




0.5 1 1,5 2 
zf solar 




-2-10 1 
[Fe/H] 



40 

h 30 

g 20 
Ed 

rt 10 





1 1 1 1 1 1 1 1 1 1 1 1 
13+5 Gyr 



^ +Skxiw 



0.5 1 1.5 2 
z/ solar 



30 



20 

a 

g 10 



1 1 1 1 1 1 1 1 1 1 1 1 

13+5 Gyr: 




-2-10 1 
[Fe/H] 



Fig. 13. Metallicity distribution in NGC 3311 assuming different age distributions. We apply the same color coding 
and method for displaying the abundance distribution as in Figure IT21 



of (Fe/H) ^-0.1, assuming a unimodal distribution of the 
metallicity. For metal-poor AND a metal rich-population 
of the same age we obtain mean metallicities of -0.27 
and +0.12 dex, respectively. For a mixed-age GCS, the 
mean metallicity of the 2 Gyr old population is derived 
as (Fe/H) 2G j, r =+0.7, with (Fe/H)i 3Gyr =-0.9 and -0.3 for 
the metal-poor and metal-rich old population. Assuming 
5 Gyr as the age of the intermediate age population, the 



results change only marginally to (Fe/H) 5( 3j, r =+0.34. The 
metallicities derived for the young metal-rich populations 
in either case are very high; compared to observed spec- 
troscopic metallicities for other second generation GCs 
(e.g. ISchweizer fc Seitzer 19981 |Goudfrooij et al. 2001al 
|Goudfrooij et al. 2001b| ) and must be regarded as very 
tentative. Several issues have to be discussed here, which 
all affect the derived metallicity distribution. 
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Fig. 14. Metallicity distribution in IC 4051 assuming a combination of a 13 and 2 Gyr old population (upper panels), 
a pure 13 Gyr old population (middle panels) and a 13 +5 Gyr old mixed population (lower panels). On the left we 
show the best bimodal Gaussian fit while on the right that for a single Gaussian. Clearly the mixed populations are 
best fitted by a double Gaussian (dashed-dotted lines). The fit to the metallicity distribution of a 13 Gyr population 
gives similarly good results for a single and a double Gaussian, although the \ 2 °f the fit is slightly better for a 
single Gaussian. In the left panel the dashed-dotted lines show the two Gaussians whereas the solid line shows their 
combination. 



— Again, the basic requisite to derive globular clus- 
ter metallicities is an accurate age determination. 
Although the combined optical and near-infrared pho- 
tometry allows us to derive the mean age of sub- 
populations, the age- metallicity degeneracy is not 
entirely broken. Even more important, for a given 
(V — ,ff)-bin the derived metallicity depends on the 
SSP isochrone. Thus applying the mean age of sup- 
populations to calculate individual metallicities causes 
additional systematic errors. 

— In our analysis we use the SSP model isochrones by 
Bruzual & Chariot IjBruzual fc Chariot 2003jl . If we 
compare the color predictions of these models with 
Vazdekis (1999) or Maraston (2001) then we find con- 
siderable differences, especially within the red popula- 
tions. Following the Bruzual & Chariot models a 5 Gyr 
old object would by ^0.1 mag redder in (V ~ I) com- 
pared to the Vazdekis models, given the same (V — H). 
This discrepancy has only little effect on the mean ages 



of globular cluster sub-populations as long as we apply 
identical SSP models for the age determination in ob- 
served and simulated samples and do not attempt to 
derive ages of individual clusters. Therefore, our pri- 
mary result of a significant number of intermediate age 
globular clusters will not change. The fact that its age 
alters with the applied SSP models between 2 and 5 
Gyrs is of no relevance since the method itself does 
not allow more precise age estimates. To quantify the 
effect of our choice of SSP models on the metallicity 
distribution we will repeat the analysis, i.e. the deter- 
mination of the cumulative age distribution and the 
modeling procedure, using the SSP model isochrones 
by Maraston (2005) (see Appendix EJ. 
— In Paper IV of this series (Hempel & Kissler-Patig 
2004) the contamination of the globular cluster sam- 
ple with unresolved background galaxies has been dis- 
cussed. Although we try to limit the contamination by 
setting limits on (V — H) and (V — /) we can not en- 
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tirely exclude it. These color limits are derived from 
the SSP model isochrones, assuming a maximum GC 
age of 15 Gyr and a maximum metallicity of 2.5 xZq 
for the younger generation of GCs. 

As already mentioned, both GCSs have been 
discussed previously (e.g. IWoodwo rth fc Harris 2 000 
|Brodie, Larsen, fc Kissle~P atig 2000|) and to further in- 
vestigate our claims a direct comparison with previous 
studies is useful. However, in the case of NGC 3311 this 
is only possible in a more general way since the study 
by Brodie, Larsen & Kissler-Patig (2000), building on the 
same optical data, includes only detections on the wide- 
field chips of HST/WFPC2, whereas combined optical and 
near-infrared data are only available for the innermost re- 
gion, covered by the PC chip. 

In IC 4051 the (V — I) distribution for objects which 
are later also detected in the infrared, is best fitted with a 
single Gaussian centered on V-I=1.17±0.005 (see Figure 
1151 left panel), which is though somewhat redder , due 
to color selection, but still in agreement with the results 
given by Woodworth & Harris (2000) (V-I)=1.106. Using 
the linear conversion by Woodworth & Harris (2000): V — 
7=0.17[Fe/H]+1.15 the mean (V — I) colors translate into 
mean metallicities of w+0.11 and w-0.3, respectively. We 
note again that these values refer to different samples, 
based on the observed field of view. 

In the discussion about the metallicity distribution 
for the NGC 3311 GCS we have to be especially cau- 
tious given the above caveats about the dependence of 
the metallicity distribution on the derived ages and is in- 
cluded rather for completeness than for solid argumenta- 
tion. Comparing the distribution of (V—I) color, as shown 
in Figure [TBI fright panel), with the results obtained for 
the wide-field chip detections, reveals as well a bimodal 
color distribution with the center of the blue and red 
sub-population at 0.91±0.01 mag and 1.11±0.01 mag, re- 
spectively. In comparison, in the work by Brodie, Larsen 
and Kissler-Patig (2000) the (V — I) color distribution 
peaks at 0.91±0.03 mag and 1.09±0.03 mag, resulting in 
a metallicity of [Fe/Hjw-1.5 and «-0.75 dex. The color- 
metallicity conversion follows the relation given in Kissler- 
Patig et al. (1998): [Fe/H] = (-4.50±0.30)+(3.27±0.32)(V- 
I). In contrast, the metallicities presented in this work, 
derived from the (V — H) color directly, but are based on 
the estimated age and a subsequent application of the SSP 
model isochrones. 

5. Summary 

We have investigated the age structure of GCs in two lu- 
minous early- type galaxies, IC 4051 and NGC 3311, each 
with large GC populations, using the combination of op- 
tical and near IR colors to help break the age-metallicity 
degeneracy. In both galaxies, we find evidence for a signif- 
icant fraction of young/intermediate age clusters («2- 5 
Gyr), in addition to their old counterparts, by comparing 
the observed color distribution with models. This result 



is confirmed by the comparison between the cumulative 
age distribution in observed and simulated systems, which 
not only indicates a second major cluster formation event 
but also allows to set some constraints on the time frame. 
Further evidence for a less passive evolution of the galaxies 
is given by independent features (dust lanes, co-rotating 
core), but especially in NGC 3311, the effects of small 
number statistics make our results tentative. Based on the 
small number of objects in this sample, our detected frac- 
tion of 50- 60% of young globular clusters in the total 
GCS in NGC 3311 has a large uncertainty. The results for 
IC 4051 are more statistically reliable, suggesting more 
than half of the globular clusters are young/intermediate 
age objects. This fraction should be regarded as an upper 
limit, considering that most of the GCs excluded from the 
age analysis are probably mainly old, and that our obser- 
vations are biased to the central regions, where younger 
objects should preferentially be found. Thus, at this stage 
the global importance of the star forming event leading 
to a second generation of globular clusters in NGC 3311 
and IC 4051 cannot be firmly established. More data, 
specially sprectroscopic, over a wider galactocentric field, 
would be useful for clarifying the nature of these popula- 
tions. Nevertheless, the improved age resolution of com- 
bined optical and near-infrared photometry reveals age 
sub-populations, which could not be detected using opti- 
cal data alone. 

Preceding the results presentet in Appendix El we 
want to include the effects of SSP model dependency 
in this summary. Although of little consequence for the 
age distribution the choice of the SSP model is becomes 
important for the derived metallicity distribution. The 
extremely high mean metallicity of the 2 Gyr old popula- 
tion ( [Fe/H]«s 0.8 dex, see Section 0J), obtained with the 
Bruzual & Chariot (2003) models, becomes much more 
reasonable when the Maraston (2005) are applied, i.e. 
[Fe/H]« 0.2. Whether this large discrepancy is based on 
the models themselves, with only the Maraston models 
including the TP-AGB phase, or is merely an artifact 
of the interpolation of the metallicity, or both, still 
needs to be investigated. Nevertheless we note that the 
difference in metallicity decreases with an inceasing age 
of the cluster population. Due to the small number of 
metal-poor objects we do not discuss model dependencies 
for this sub-sample here. 



Acknowledgments 

Support for this work was provided by NASA through 
grant number GO-07280 from the Space Telescope 
Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under 
NASA contract NAS5-26555. K.Mighell kindly provided 
advice concerning the CCDCAP package. D.G. would like 
to thank Nick Suntzeff and Malcolm Smith of the staff of 
Cerro Tololo Inter- American Observatory for their kind 
support of this project. D.G. gratefully acknowledges 



Hempel et al.: Extragalactic Globular Clusters in the Near Infrared V 



17 






25 - 




20 : 


= CD 


1 5 : 


= _Q 


I E 




- ^ 
= z 


10 : 




5 : 




: 



0.6 0.8 1.0 1.2 
V— I [mag] 



.4 1 . 




0.4 0.6 0.8 1.0 1.2 1.4 1. 
V— I [mag] 



Fig. 15. (V - I) color distribution for globular clusters in IC 4051 (left) and NGC 3311 (right) detected on the PC 
chip, for which as well H-band data are available. The open histogram represents the complete VIH sample, whereas 
the shaded histogram includes only objects for which the photometric errors, 6(V — I) and S(V — H) have been derived 
to be <0.15 mag. The solid curve shows either the single Gaussian (for IC 4051) or the sum of the two Gaussian 
functions (dotted and dashed lines) in case of NGC 3311. 



support from the Chilean Centro de AstroHsica FONDAP 
No. 15010003. The authors wish also to thank the 
anonymous referee for her/his comments, which helped 
immensely to improve the manuscript. 



Appendix A: SSP model dependency 

In the recent past, tremendous progress has been made 
in the understanding of stellar evolution and one of the 
major outcomes is more and more accurate models of how 
colors and spectral features for single stellar populations 
evolve in time, depending on metallicity and other pa- 
rameters, e.g. the structure of the horizontal branch (e.g. 
IMaraston &: Thomas 20"00l |Lee, Lee fc Gibson 2002"| ). 
Although this paper does not intend to argue in favor of a 
specific SSP model wc would like to show how the results 
depend on our choice of the SSP model. In the main part 
of this paper we apply the new version of the Bruzual & 
Chariot model isochrones (2003). As shown in Hempel 
& Kissler-Patig (2004) and Hempel (2004) the choice 
between the Bruzual & Chariot models (2000, 2003) or, 
alternatively, the model by Vazdekis (1999) or Maraston 
(2001) affect the results only marginally, i.e. the detection 
of intermediate age clusters is essentially independent 
of the selected model. However, in order to investigate 
the subsequent advances made in SSP modeling we 
will as well conduct our analysis applying the latest 
version of the SSP models by Maraston (2005), assuming 
a red/intermediate horizontal branch structure and a 
Salpeter IMF. The (V - I) and (V - H) colors for both 
the Bruzual & Chariot (2003) and Maraston (2005) model 
are compared in Figure IA.1I We note that unlike earlier 
releases of both SSP models the color predictions for a 
given age do not differ significantly, except for metal-poor 
objects with ages /leq2 Gyr. Prior to our analysis we 
compared as well the model predictions for the (V — I) 
and (V — H) color assuming a red/intermediate and 
a blue horizontal branch and various ages (> lGyr) as 
offered by the SSP models. With the exception of a 1 



Gyr old stellar population, the change in color is smaller 
than the photometric error in the observational data 
and in fact smaller than the color difference between 
models using different IMF's, especially within the red 
(V — H) color range where our technique is sensitive to 
the detection of young /intermediate GCs. 

The color-color diagrams, for both globular cluster sys- 
tems are again shown in Figures IA.2I and I A. 31 In contrast 
to Figure [3] and [7| the colors are now compared to the SSP 
models by Maraston (2005). Still the color distributions 
show a large fraction of globular clusters with a (V — I) 
color which is too blue to be consistent with a population 
older than 5 Gyr. 

In order to detect GC sub-populations within a 
galaxy we derive the cumulative age distribution (see 
|Hempel et al. 2003| |Hempel fc Kissler-Patig 2004| | based 
on optical and near-infrared colors in the same way as in 
Section 3.2. Hereby the GC density with respect to the 
model isochrones is of more importance than the color of 
each GC separately. As described before, the GC density 
is determined as the number of objects with a (V — I) 
color redder and therefore older than a given isochrone. 
Previously those number counts were estimated for the 
1,2,3,5,7,10,13 and 15 Gyr isochrone. Using the new 
Maraston models requires a small alteration in this 
procedure. Despite the non-monotonic pattern of the 1 
Gyr isochrone we follow the same counting procedure, 
taking advantage of the fact that, independent from 
metallicity, the model predictions for (V — I) give a 
redder color for a 2 Gyr old population compared to 1 
Gyr population. Nevertheless, in the simulations we as- 
sume the second generation of GC to be at least 2 Gyr old. 

Comparing the cumulative age distributions in both 
observed systems with our set of simulations, shows at first 
significant differences, if the Maraston SSP models are ap- 
plied (Figures I A . 41 and I A . 5|) . Hence we need to emphasize 
again that this approach does not aim at high precision 
age estimates but is rather a tool to detect significant sub- 
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Fig. A.l. Comparison between the color predictions of the Bruzual & Chariot (2003, dotted line) and the Maraston 
(2005, solid line) SSP models. Clearly the largest discrepancy is found for objects younger than 3 Gyr, whereas for 
older populations the difference in (V — I) is in the range of the photometric error or below. We note that the models 
are given for different metallicities (Bruzual & Chariot: -2.25,-1.35,-0.64, -0.33,0.0,-1-0.56 dex; Maraston:-2. 25,-1. 35,- 
0.33,0.0, +0.35 dex). 
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Fig. A. 2. (V — /) vs. (V — H) color-color diagram for IC 4051. As in FigureGlonly GCs with photometric errors <0.15 
mag have been selected. The 1,5, and 15 Gyr isochrones are marked as solid, dotted and dashed line, respectively. Here 
the SSP model by Maraston has been applied. 



populations of GCs. Ideally, a successful detection should 
be complemented by spectroscopic investigations to de- 
termine the age more accurately. Due to the large dis- 
tance of most early-type galaxies (especially in high den- 
sity environment) spectroscopic surveys are hardly feasi- 
ble. Keeping this in mind we find that the conclusions 
about a second population of intermediate age globular 
clusters with the selected sample is confirmed, although 
the relative size of the second GC generation is shifted to 
larger values, i.e. using the Maraston 2005 models leads 
to an even larger fraction of intermediate age clusters. 
Since the latter is affected by various biases (color cuts, 
spatial bias, etc.) we consider large numbers of intermedi- 



ate/young GCs to be an upper limit at best. Considering 
the importance of horizontal branch structure for the mod- 
eling of SSP isochrones this result seems to be surprising. 
On the other side we have to be aware that a major impact 
of the HB structure is expected for the colors of metal- 
poor objects (see Figur( lA~T|l . which we exclude from our 
analysis by setting color cuts. 

To round things up we will finish this Section by 
deriving the metallicity distribution in IC 4051, following 
the procedure described in Section but now based on 
the Maraston (2005) models. As before we calculate the 
metallicity under the assumption of the globular clusters 
being 2, 5, or 13 Gyr old, depending on their (V — I) 
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Fig. A. 3. (V — I) vs. (V — H) color-color diagram for NGC 3311. As in Figure |3] only GCs with photometric errors 
<0.15 mag have been selected. The 1,5, and 15 Gyr isochrones are marked as solid, dotted and dashed line, respectively. 
Here the SSP model by Maraston has been applied. 
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Fig. A. 4. x 2 -test result for IC 4051. Different levels represent the x 2 -values of the comparison between the cumulative 
age distribution in the selected GC sample and a set of 66 simulated systems. For both, the modelling and the 
determination of the age distribution the SSP models by Maraston (2005) have been applied. Similar to the application 
of the Bruzual & Chariot models (2000) the best fitting models contain a large fraction of young/intermediate age 
globular clusters. 



color with respect to the corresponding isochrones. The 
metallicity [Fe/H] for each objects is then interpolated 
for a given (V — H). 

Comparing the color predictions of both models we 
note that stellar populations older than 3 Gyr are al- 
most unaffected, the Bruzual & Chariot (V — I) colors 
being somewhat bluer than their Maraston counterparts. 
Translated into ages the globular clusters are slightly 
younger when the Maraston models are applied. There is, 
however, a significant difference in (V — H) with respect to 
the corresponding metallicity values. As shown in Figure 
IA. li the (V — H) color of a 1 or 2 Gyr old population of low 



metallicity (i.e. [Fe/H]«-2.3) is much bluer in the Bruzual 
& Chariot models. Consequently, the metallicity distri- 
bution for the IC 4051 GCs will be shifted toward lower 
metallicities when the Maraston models are applied. This 
can be seen in the left panels of Figure lA~6l The metallic- 
ity for 2 Gyr old objects (red histogram) peaks at much 
lower values compared to Figure El left panels. For the 
first generation of GCs (13 Gyr) the metallicity estimates 
change only marginally and are also similar for a 5 Gyr 
population within the anticipated accuracy of 0.2 dex. 

Based on the above results we therefore conclude that 
our primary goal- the detection of age sub-populations 
in GCSs- is not significantly influenced by the choice of 
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Fig. A. 5. Contour plot representing the results of the x 2_ test comparing the cumulative age distribution in the 
NGC 3311 GCS (selected sample) and the Monte-Carlo simulations. As in Figure lA~4l the SSP models by Maraston 
(2005) have been used. 
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panels). The different panels and the color coding of the histograms correspond to Figure IT51 



the SSP model. Nevertheless, we note that our method is 
somewhat handicapped when the SSP model isochrones 
are not monotonic, e.g. for a given age an unique cor- 
relation between the two colors can not be established. 
There is, however, no doubt that state of the art SSP 
models should be applied, especially when the informa- 
tion we want to extract from the globular cluster colors 
become more and more detailed (including metallicity dis- 
tributions) and an increasing accuracy of the data allows 
more precise calculations. 
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